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We studied the hadronic decays of excited states of heavy mesons {D, Ds, B and Bs) to lighter states by 
emission of tt, t] or K. Wavefunctions and energy levels of these excited states are determined using a Dirac 
equation for the light quark in the potential generated by the heavy quark (including first order corrections in the 
heavy quark expansion). Transition amplitudes are computed in the context of the Heavy Chiral Quark Model. 



1. INTRODUCTION 

In the context of the most general Quark 
Model, a heavy-light meson {H) is modeled with a 
light quark [q) bound to a static source of chromo- 
electro-magnetic field (the heavy quark h). Its 
hadronic transitions can be computed assuming 
that only the light quark enters in the reaction 
through and effective coupling qa (Heavy Chiral 
Quark Model [|). 

Our work generally follows that of ref. 
We differ in the choice of the potential. Moreover, 
we included for the first time the mixing effects 
in the spectrum and decay amplitudes. 

2. SPECTRUM 

2.1. Notation 

The Dirac wavefunction, $, of the light quark 
can be determined by solving the eigenvalue prob- 
lem 
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where H is the Hamiltonian of the system and 
^ is 4-spinor that represents the wavefunction of 
the system. In our notation 
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where n is the radial quantum number, (, is the 
orbital quantum number, j, m are the total spin 



of the light quark and its z component, J, M are 
the total spin of the meson and its z component, 
/",/^ are radial wavefunctions and Y^ are the 
usual spherical harmonics. ^5 is the 2-spinor as- 
sociated to the heavy quark h and S is its spin. 

Our convention for the phase and the normal- 
ization is such that 
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*TaIk presented by 



2.2. Choice of the potential 

Ignoring l/ruh corrections, the most general 
form of the Hamiltonian that appear in eq. (|l|) 
is 

where Vs is a spin independent potential, Vy is a 
spin dependent potential and Mh is a total en- 
ergy shift (not to be confused with the mass of 
the heavy quark mh that appear in corrections to 

Asymptotic freedom suggests that at short dis- 
tances Vy ~ 1/r dominates, while lattice simula- 
tions suggest that at large distances Vs ^ r domi- 
nates. As it was observed in ref. ||], the choice of 
Coulomb-like potential at short distance is incon- 
sistent with l/rrih spin-dependent correction to 
iJ", because of ultraviolet divergences. The solu- 
tion of the problem is that, in the context of the 
Dirac equation with a finite mh-, it is not correct 
to localize the heavy quark with a delta function 
since one must take into account the spatial de- 
grees of freedom of the heavy quark. Our prag- 
matic approach to the problem is that of delocal- 
izing the heavy quark within a length scale 1/A 



assuming a Gaussian wave-function, $(a;), for the 
former. The effective potential felt by the light 
quark is, therefore, a convolution of the Coulomb- 
like potential with the square of the wave- function 
of the quark: 
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Our choice for the spin-independent part of the 
potential is 



Vs (r) — br A- c 



(6) 



(notice from eq. (0) that c is not a physical param- 
eter since it can be re-absorbcd into the definition 
of niq). 

2.3. 1/m/i correction 

We solve the eigenvalue problem associated to 
the eq. (|l|) using the Hamiltonian in eq. (Q) and 
the potentials (||,||) . In this way we determine the 
radial wave- functions /^^, /^^ and the associated 
eigenvalues Entj- 

1/mfi corrections to the the Hamiltonian have 
been derived in [Q in the Bethe-Salpeter formal- 
ism. They are responsible for the spin-orbit inter- 
action, the hyperfine spitting and the mixing of 
states with the same j. We include these effects as 
a perturbation to the energy levels (SEnijj) and 
also determine the mixing coefficients for each 
doublet of states. 

In terms of the /' the l/rrih correction to the 
energy levels reads as 
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where one can identify the contributions propor- 
tional to p^: 
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the spin orbit interaction, for j = i - 
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Figure 1. Spectrum of excited B mesons. 
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Table 1. Tabulated spectrum for the observed 
states used in the fit together with predictions 
for some excited B states. All units are in MeV. 
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Figure 2. Orbitals for some excited B mesons. 



and the hyperfine splitting contribution: 
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2.4. Choice of parameters and predictions 

The nine parameters of the model (a^, A, 
b, mq\q=u, ruqlq^s, rnh\h=c, Mh\h=c, rnh\h=b-, 
Mh\h=b) are determined by best fit (minimum x^) 
to the known spectrum of excited states of the D, 
D,, B and B^ mesons. Our results are: 



as = 0.339 






A == 2.820 GeV 






b ^ 0.257 GeV 






m„ ^ 0.073 GeV 


ms 


= 0.214 GeV 


ruc = 1.52 GeV 


Me 


= 1.52 GeV 


mb = 4.67 GeV 


Mb 


= 4.68 GeV 



The spectrum for some of the states plotted in 
fig. ^ and tabulated in table 1. A density plot of 
\fltjYa? and \fl,^jY^'~^? for some excited states 
of a B meson is reported in fig. El 

2.5. Tests of the model 

Despite the good fit of the mass spectrum we 
decided to test our model by comparing some of 



the transition amplitudes with some recent lattice 
results. In particular we computed 
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using our chiral quark model (where the only un- 
known parameter is the overall normalization, gA, 
which is the effective coupling of the quark to the 
axial current, A^) and comparing it with the lat- 
tice result of ref. Q . The comparison is shown in 
fig. ||. In the plot the point at r = is used 
to fix the relative normalization. A more so- 
phisticated analysis of the lattice results gives^ 
5A == 0.42 ± 0.09. 

3. HADRONIC DECAYS 

3.1. Decay amplitudes 

We consider here the most general hadronic 
transition 



H' -^ H 



(13) 



where H' is an heavy meson with associated wave- 
function 4"', iJ is an heavy meson with associated 
wave-function \E' and a: is a light meson with mo- 
mentum p. Such a transition is mediated by a 
matrix element of the form 
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where X is the 4x4 spin matrix that characterize 
the transition. In the particular case in which 
the light meson x is a pseudoscalar (tt, ry or K) 
X = ^^^, while if the light meson x is a vector (p, 
uj or K*) X — ^ (and e is the polarization vector 
of the outcoming vector meson) . 

The exponential can be expanded in products 
of spherical harmonics and spherical Bessel func- 
tions (j'fc), thus giving 
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We computed Ai^ {X, p) for a complete set of spin 
matrices X and proved that it can alway be re- 
duced to integrals of the form 
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where cl (X) depend also on the quantum num- 
bers of the mother and the daughter mesons but 
not on the radial wavefunctions. 

3.2. Decay widths 

The decay width for transition of eq. (O) , when 
light meson x is emitted in an eigenstatc of the 
total momentum p and of the angular momentum 
4, is given by 
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where C = VS, l/'s/Sj 2/V3 or 1 for tt, tj (for a 
nonstrange heavy meson), -q (for a strange heavy 
meson) or K respectively; M', J' and the mass 
and total angular momentum of the mother me- 
son H'; AI, J are the mass and angular mo- 
mentum of the daughter meson H, and p can 
be determined by energy-momentum conserva- 
tion (po = M' - M,p=\p\ = Vp^~^). The 
total decay width and the branching ratios are 
defined as 
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where the sum on x spans all the hadronic decay 
modes with emission of a light pseudoscalar me- 
son. The total width for some of the mesons is 
reported in table 1. 

4. CONCLUSIONS 

We computed the spectrum and the width of 
hadronic decays of excited D, Dg, B and Bg 
mesons in the context of the chiral quark model. 
As an example, we report here some hadronic de- 
cay channels for the first radial excited B meson 
(not including the p decays) 
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Br = 77%) 
Br = 22%) 
Br ^ 100%) 
Br = 0.33%) 
Br = 0.03%) 



More complete tables will be published on a 
separate paper. 




Figure 3. Comparison between a prediction of our 
model and the lattice QCD result B. 



We finally remark how our model is able to 
fit the masses of observed excited states within 
less than 10 MeV discrepancy (within 4 MeV in 
average) better than was done in preceding works. 

This work was performed at Fermilab, a U.S. 
Department of Energy Lab (operated by the Uni- 
versity Research Association, Inc.), under con- 
tract DE-AC02-76CHO3000. 
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